Observation of Sommerfeld precursors on a fluid surface 
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We report the observation of two types of Sommerfeld precursors (or forerunners) on the surface 
of a layer of mercury. When the fluid depth increases, we observe a transition between these 
two precursor surface waves in good agreement with the predictions of asymptotic analysis. At 
depths thin enough compared to the capillary length, high frequency precursors propagate ahead 
of the "main signal" and their period and amplitude, measured at a fixed point, increase in time. 
For larger depths, low frequency "precursors" follow the main signal with decreasing period and 
amplitude. These behaviors are understood in the framework of the analysis first introduced for 
linear transient electromagnetic waves in a dielectric medium by Sommerfeld and Brillouin jlj. 
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One feature of linear wave propagation in a dispersive 
medium is the existence of precursors (or forerunners). 
This terminology traces back to the fact that they gen- 
erally arrive sooner than the "main" signal. This tran- 
sient response is due to the propagation of the fastest 
high frequency components of the initial spectrum. Al- 
though predicted as early as 1914 by Sommerfeld and 
Brillouin P], experimental observations of forerunners 
are very few and only qualitative, mainly dealing with 
electromagnetic (e.m.) waves in dielectric medium in 
the microwave or optical Q frequency range. Such 
Sommerfeld forerunners have also been predicted in vari- 
ous dispersive media such biological 0| or viscoelastic 
ones, and have been recently shown to be linked to the 
non violation of Einstein causality during superluminal 
light pulse propagation in the region of anomalous disper- 
sion 6J. However, for waves in fluids, the observation of 
forerunners is still lacking despite some effort performed 
with acoustic waves in superfluid '^He 7j or with pressure 
waves in fluid-filled collapsible tubes . We report here 
the first observation of two types of Sommerfeld forerun- 
ners, which can coexist or not, on a thin layer of mercury. 
The non-monotonous dispersion relation of waves on the 
surface of a fluid lead to a rich variety of such transient 
wave phenomena. In the shortwavelength limit, when 
capillary effects are dominant (the analogue to "anoma- 
lous" dispersion jof for e.m. waves), only the fast high fre- 
quency components of the initial excitation are observed 
arriving before the main pulse, the so-called Sommerfeld 
precursor for e.m. waves. This transient is characterized 
by small amplitude and rapid oscillations (with respect 
to the main signal). For longer wavelengths, when grav- 
ity is no longer negligible ( "normal" dispersion) , a slower 
low frequency precursor is also observed (no analogue 
exists for e.m. waves). The forerunners are found in 
good agreement with the predictions of asymptotic anal- 
ysis based on the "stationary phase method". We note 



that this study also allows us to connect the forerunner 
concept of electromagnetic waves to the well-known hy- 
drodynamic transient surface wave phenomena |lOl |. and 
their applications to submarine eruptions [1^ . 

The experimental setup consists of a 1.5 m long hor- 
izontal PMMA channel, 7 cm wide, filled with mercury 
up to a height, h, where: 2 < /i < 14 mm. h is mea- 
sured to a precision of ± 0.02 mm by means of a depth 
gauge using a micrometric linear positioner. The proper- 
ties of the fluid are, density, p = 13.5 lO'^ kg/m"^, dynamic 
viscosity, rj — 1.5 10~^ Ns/m^ ^^'^ surface tension 
7 — 0.4 N/m |l3j . Surface waves are generated by an im- 
pulsional excitation provided by the horizontal motion of 
a rectangular plunging PTFE wavemaker driven by an 
electromagnetic vibration exciter. They are generated 
10 mm inward from one end of the channel and the local 
displacement of the fluid in response to this excitation is 
measured simultaneously by a nonintrusive inductive sen- 
sor and by an optical technique The inductive sen- 
sors, 3 mm in diameter, are suspended perpendicularly 
to the fluid surface at rest. The linear sensing range of 
the sensors allows distance measurements from the sensor 
head to the fluid surface up to 2.5 mm with a 5 V/mm 
sensitivity. An optical determination of the local slope of 
the surface is also performed. Using a position sensitive 
detector, we have recorded the deflection of a laser beam 
by the surface wave; the computation of the surface el- 
evation from the optical signal is in excellent agreement 
with the direct inductive measurement of the shape of the 
wave [see inset of Fig.lSJb-c)]. Both sensors are mounted 
on a horizontal linear positioner at a distance x from the 
wavemaker, < x < 1.2 m. The optical and inductive 
methods are complementary as the spatial resolution and 
the sensitivity of the optical technique are higher, but the 
inductive technique provides a direct measurement of the 
surface displacement and does not require signal process- 
ing. Both techniques are not limited by their response 
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time in the frequency range of surface waves. The choice 
of mercury has been motivated by the possible use of the 
inductive measurement technique and also because of its 
low kinematic viscosity which is an order of magnitude 
smaller than that of water, thus strongly reducing wave 
dissipation. Photographs of the free-surface of mercury 
are made by means of a camera mounted above the center 
of the channel. A typical pattern is displayed in Fig. 
showing the high frequency Sommerfeld precursors, Sh, 
ahead of the main signal followed by low frequency ones 
Sl- 




'H 



FIG. 1: Photograph of typical wave precursors at the surface 
of mercury (seen from above). Pulse fronts are located on 
the left. Fluid depth is h = 3.7 mm (Bo=0.22), and the full 
vertical scale corresponds to the 7 cm canal width. 

To understand this behavior, one can introduce the 
dispersion relation for surface waves, neglecting dissipa- 
tion. 




tanh kh 



(1) 



with Lo the wave pulsation, k the wave number and g the 
acceleration of gravity. From Eq.J^I, we can define the 
capillary length, Ic = ^/t/ (pg), and the Bond number. 
Bo = {Ic/hy^. In the long wave length approximation 
or "shallow water" limit [kh << 1), Eq. ^ may be ex- 
panded, and the group velocity becomes. 



dk 
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gh[l ~ a2k h 
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with 02 = (i - Bo), and 04 = (19/90 - Bo/2 - Bo'^/3). 
Equation ^ shows that a minimum of the group velocity 
exists only when < Bo < 1/3, Bo = 1/3 corresponding 
to a critical depth, he- FigureElshows a qualitative sketch 
of the wave group velocity as a function of fc, for Bo above 
and below 1/3. As explained hereafter, the existence of 
this minimum has a strong influence on the dispersion of 
an initial disturbance. 

Generally speaking, if we look at the propagation 
along Ox of an initial perturbation Co{x) in a disper- 
sive medium, ^(a;, t) (e.g. the free-surface deforma- 
tion) is formally given by the Fourier Integral ({x, t) = 
I-x C,o{k)e^'^^dk, where = kx/t — ui{k), with uj{k) solu- 
tion of Eq. and where Co{k) is the Fourier transform 
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FIG. 2: Group velocity vs wave number for gravito-capillary 
waves. At fixed x and t, the stationnary phase method 
[Vg{ks) = x/t] leads to solutions with wave number kg. For 
h > he (Bo < 1/3): at small t, high frequency Sommerfeld 
precursors, Sh, are predicted ahead of the "main signal" with 
a period which increases with time (fesj, decreases). When t 
is large enough {x/t < \fgh) Sh switches to a low frequency 
Sommerfeld precursor, Sl, with a decreasing period (ks^ in- 
creases). A Brillouin precursor corresponds to a constant pe- 
riod {kmin)- The typical signal amplitude vs time is sketched 
on the right. The case h < he (Bo > 1/3), where only the 
fast Sh solution exists, is displayed in the inset. 



of C,o{x) The method of the stationary phase [l^ll5| 
is particulary useful for the asymptotic behavior of these 
integrals: for both large x and t {x ^ L ~ 2 cm and 
t':S>T~ 100 ms, the size and the duration of the initial 
disturbance), with x/t held fixed (for an observer trav- 
eling at this given speed), the main contribution to the 
integral is from the neighborhood of stationary points kg 
such that 



dk 



,i.e. 



dhJ 

dk 



= Vg{ks) = x/t 



(3) 



the other components oscillate too rapidly in order to 
contribute. One can apply graphically this stationary 
phase method to the dispersion relation of Eq. (QJ, as 
shown in Fig. [2 At a fixed point of observation x, the 
main contribution to the surface deformation ^(x, t) at 
any time t results from the points on the group velocity 
curve equal to x/t (see Fig. For h > he (Bo < 1/3), 
three types of precursor are predicted: the fastest signal 
is the high frequency Sommerfeld precursor, Sh, (from 
the capillary branch) ahead of the "main signal" (arriving 
at to with velocity x/to = y/gh); then, the low frequency 
Sommerfeld precursor, Sl (from the gravity branch), and 
finally, at , the so-called Brillouin precursor (minimum 
of the Vg{k) curve, i.e. (p" = 0) as defined in the frame- 
work of e.m. waves. For h < he (Bo > 1/3), only the 
fast Sh solution exists (see Fig.EJ. Fig.[21also shows that 
Sh has an increasing period as time goes on (contrary to 
Sl), whereas the Brillouin precursor has a constant pe- 
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riod. This rich variety of such transient waves comes 
from the non- monotonous dispersion relation of Eq. 
Finally, returning to the general case, one can have ac- 
cess to the disturbance profile C{x,t) which is given by 
the above Fourier Integral. A Taylor serie expansion of 
(j){k) in the neighborhood of kg leads to 0, Hsf 



2tt 



dvg/dk\. 



cos [(l>{ks)t ± 7r/4] . 



(4) 

Thus, the shape of the precursor notably depends on 
the Fourier transform of the initial disturbance, (^o{ks), 
whereas the temporal evolution of the precursor period 
depends only of the dispersion relation through (j){ks). 
Therefore, the measured period can be easily compared 
with the theoretical prediction, whereas the prediction of 
the experimental profile requires the precise knowledge of 
the initial conditions. Note that, at the minimum of the 
group velocity where the Brillouin precursor is predicted, 
Eq. 0) is no longer valid, and the correct asymptotic be- 
havior is found by keeping higher orders in the Taylor 
series for ^(fc) ^E^. 

We have performed an experimental study of precur- 
sor waves for a fluid layer height in the range 2.12 < 
h < 13.75 mm; thus 0.02 < Bo < 0.67. For mercury 
{Ic — 1.74 mm), the critical case Bo = 1/3 corresponds 
to /ic ~ 3 mm. A horizontal impulsion is imposed to 
initiate the surface wave. In response to this initial dis- 
turbance, the free-surface profile is recorded at a fixed 
distance from the wave generator, and displayed in Fig. 
Ela) for h < he (Bo > 1/3) in Fig. Olb) for h > he (Bo 
< 1/3). Figure |2Ia) shows the high frequency Sommer- 
feld precursor which propagates at least two times faster 
than the main signal (t/to > 0.5), with an increasing 
period (defined as the time between two successive ex- 
trema) as graphically predicted in the inset of Fig. |21 
For a deep enough fluid (Bo < 1/3) this fast precursor 
coexists, as shown in Fig.OIb), with low frequency Som- 
merfeld precursor which appears behind the main signal, 
with a decreasing period as graphically predicted in Fig. 
121 Moreover, the amplitude of the fast precursor, Sh, 
is much smaller than Sl (see vertical scales in Fig. 01 a- 
b)). This can be understood by assuming an initial nor- 



malized hump of section s, Co{x) = /{s^ + x^), whose 
Fourier transform is Co(fc) ~ e~^^ . e^'"' being maximum 
for fc = 0, the precursor Sl with the larger spatial scale 
has the larger amplitude. For the same depth of fluid as 
in Fig.Ofb), Fig.Ofc) shows the profile of waves after 0.6 
m of propagation. The fastest precursor has disappeared, 
and only the contribution due to the gravity branch, i.e. 
Sl, is observed. This is linked to the arguments given 
above, and also to viscous dissipation. The inset of Fig. 
inib-c) shows good agreement between optical and induc- 
tive measurements, except near the front wave where the 
small and fast Sh forerunners are not resolved by the 
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FIG. 3: Free-surface profiles of Sommerfeld precursors as a 
function of dimensionless time (to = xj^fgli): (a) Type Sh for 
h = 2.12 mm at a; = 0.2 m, (b) Type Sh and Sl for h = 7.2 
mm at j: = 0.2 m, and (c) Type Sl for h = 7.2 mm, x = 0.6 
m. Inductive (thick line) or optical (thin line) measurements. 
Pulse fronts are located on the left. Insets show a comparison 
between both techniques. 
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inductive method (see inset of Fig. Ufb)). 

Finally, the period of each forerunner is measured all 
along its propagation (0.2 < x < 0.8 m) by recording 
the time between successive maxima of the amplitude. 
The periods are displayed in Fig. 0] as a function of to /t 
(to = x/\/gh) for various fluid depths h corresponding 
to 0.02 < Bo < 0.67. For each height corresponding to 
Bo > 1/3, the Sh precursor period increases as the waves 
propagate. For < Bo < 1/3, both Sh and Sl precur- 
sors are observed, and the period of Sl decreases with 
time. We can map Fig. 21 on Fig. [3 with a 90° rotation 
and {x/t)/y/gh can then roughly be viewed as the group 
velocity curve of the stationary mode ks as a function of 
l/ks- Therefore, Fig. 0] shows that Sh precursor veloci- 
ties are supersonic {vg{ksH) = x/t > \/gh), whereas Sl 
precursor is subsonic (i'g(fcs^) < \/gh). For each h, all 
the data recorded at different x, lie on a single curve pre- 
dicted by Eq. (Q), which is the parametric plot of 27r/a;(fc) 
as a function of Vg{k) = du/dk for various values of k. 
Note that a Brillouin forerunner is never observed in our 
experiments, and Sh disappear for large h (absence of * 
and o-marks on Fig.0J. Note also the absence of Sh for 
{x/t)/y/gh < 1. Although they are predicted from Fig. 
El they are much smaller than Sl, and thus cannot be 
observed when they travel at the same velocity. 



300- 



250 

-.200 
in 

|l50 

'03 
Cl 



100 

so 





Lcufii frequency / / r / 
Sonitiierf&lcl's pcecuiMts/ */ / / 


1 i 

1 High (reqiiencv 


/ /" hi 




h 1 1° 




/ h 1^ 
1*1 i \ 








r r r k 




Y i £^> f^ 
1 f> y 




- - - ^ L nip rftLJ ° 




BriHoLiln -L^^-J' "p l 
pr-eojFsors V-t — ^"H 






1^ - 



0.5 



1 

(x/t)/(qll)^' 



1.5 



FIG. 4: Period of Sommerfeld precursors Sh and Sl as a 
function of {x/t)/^/gh for various heights ft = 2.12 for depres- 
sion (v) or elevation (A) pulses, 3.4 (x), 5.6 (□), 7.2(o), 10.4 
(o), 13.75 (*) mm with 0.2 < a; < 0.8 m. For each value of h, 

the theoretical Sommerfeld (— ) and Brillouin ( ) precursor 

periods are extracted from Eq. Q (see text for details). 

In conclusion, we have reported the observation of 
Sommerfeld forerunners in the shallow water limit of sur- 
face waves and found that their period is in good agree- 
ment with theoretical predictions. An extension of this 
work, much easier to study in the context of hydrodynam- 
ics than in the one of electromagnetism, is to understand 
how the dynamics of precursors are changed when the 



main signal amplitude is increased such that nonlinear 
effects become important. 
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